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^ ! ABSTRACT 
00 

Context. The lag-luminosity relation (LLR) provides a way of estimating GRB luminosity by measuring the spectral lags between 
. different energy bands. 

^Nj Aims. We want to understand the origin of the LLR and test its validity. This appears especially important if the LLR is to be 
, used as a distance indicator. 

C ■ Methods. We perform a linear analysis of the lag between two spectral bands. The lag is obtained as the time interval between 
' the maxima of a given pulse in the two bands. 
Results. We get a simple expression for the lag, which shows in a very simple way how it is related to the spectral evolution of 
rS i' the burst via the variation of the peak energy and spectral indices. When this expression is coupled to the Amati relation, it 
leads to a LLR that agrees with the observational results only if the burst's spectral evolution is limited to a decrease in peak 
energy during pulse decay. However, when the variation of the spectral indices is also taken into account, the predicted LLR 
differs from the observed one. 

Conclusions. We briefly discuss some ways to solve this problem, such as a possible correlation between pulse spikiness and burst 
• • luminosity. 
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1. Introduction lag Ati^ anticorrelates with burst luminosity and propose 

the following power law relation 

L = L310^3 (Aiig/O.Ol s)"i-i^ erg.s-i . (1) 

The problem of the distance to gamma-ray bursts re- r^^^ ^^.^^ LLR was then investigated by Kocevski 

mained unsolved until the discovery of the afterglows ^iang (2003), Ryde (2005), and Ryde et al (2005), 

by Beppo-SAX. Redshifts are now obtained from optical i, r j u j i ^ j-r. 

r , r 1 , , who tound that the observed lags are a consequence 01 the 

spectra or the afterglow itself or oi the host galaxy when i ; , ^ i i j-- t ^u- j. -i, j.- c 

f , , , r burst s spectral evolution. In this contribution we perform 

the afterglow has faded away. Using the known redshifts, ^ ^^^^^^.^ ^^.^^ 1^^^ ^.^^^^^ 

it became possible to calibrate relations by linking abso- ^lR may be linked to the Amati relation and possibly 

lute burst outputs (luminosity or total radiated energy) ^^^^ variability-luminosity relation, 

and quantities directly available from the observations in 
gamma-rays. A Cepheid-like relation between variability 

and luminosity was proposed, for example, by Reichart et 2. Count rates in different energy bands 

al (2001). More recently, Atteia (2003) used the Amati re- „^ , , , rr. i i 

lation (Amati et al 2002) to introduce "pseudo-redshifts" ''"T ^ ^P^^*^^^ ^f^^ f^*^ assume a spec- 

1 . , ,j 1 r- 1 . -11 -J j-r 1 ■ 1 /-irm trum shape consisting of two smoothly-connected power 

which could be useful to rapidly identity high-z GRBs , „^ . 

r . i- 1 T j-T,- laws of respective slopes a and p at low and high energy 

from their gamma-ray properties alone, in this paper we f / \ . i i r i 

, , .1 1 1 • 1 /TTT^^ J- (Band et al, 1993). The count rate 7V,i(t) in band LB,-, i?i 

concentrate on the time lag-lummosity relation (LLR) dis- ^ ' ' ■'^ ' l o 

covered by Norris et al (2000). The lags were computed by 

Norris et al using the burst profiles in BATSE band 1 (20 

- 50 keV) and 3 (100 - 300 keV). They find that the time 



{t)=A{t) f ' B^(3{x)dx (2) 

J Xi 
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where the function A{t) depends on time alone. The li- 
mits of the integral are Xij = Ei j{l + z)/Ep{t), z being 
the redshift of the source, Ep{t) the peak energy of the 
instantaneous spectrum (in source rest frame) and Bap{x) 
the spectrum shape. 

Considering now another spectral band [Ej^, Ei], we 
can relate Nki [t) to A^^ {t) in the following way 



Nki{t) = Nij{t) 



Q Bcp{x)dx 



= Nij{t) X Tijki [Ep{t),a{t),p{t)] 



(3) 



where J^ijki can be seen as the "spectral correction" be- 
tween bands [Ei, Ej] and [Ek, Ei]. We simplify the nota- 
tion by considering only BATSE bands 1 [20, 50 keV] and 
3 [100, 300 keV] so that we have 



7V3(t)=iVi(0x.Fi3 {Ep,a,(3] 
with 

p300(l+2)/£;p 



:Fi3 = 



JlOO(l+2)/Bp '^ap\X)UX 



.50(l+.)/Bp „ . . 
J20(l+z)/Bp l^aliyxjad, 



(4) 



(5) 



We then assume that a given pulse in the burst profile 
reaches its maximum at a time ti (resp. ts) in band 1 
(resp. 3) and we estimate the lag by the difference 



(6) 



Since in most cases the observed lags are small compared 
to the pulse duration, we evaluate Atis from a linear 



(7) 



analysis of the pulse shape around ti. Using 
Ni{ti) =0, we can write 



while the spectral correction gives to same order 

^13{t) = J^i3(ii)+:^i3(ti) {t-tl) + ^fil3{ti) {t-ttf , (8) 

with .^13(^1) and ^-"13(^1) being related to the partial 
derivatives of ^13 with respect to Ep, a, and For Tis{t{) 
we have 



:^i3(ii) 
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dE^ 



ti da 



a{ti 



that they are essentially negligible for the final numeri- 
cal results. Solving Eq.(lO) to get t^, such as Ns^t^) = 0, 
finally yields 



13 



df3 



while .^^13(^1) contains nine terms. We now compute the 
logarithmic derivative of N3{t) to the first order in (t — ti) 



N^it) mit) J'lait) Niih) 
^13 (ti) _ / -^isfto V 



(t-h). (10) 



The bracket contains many terms involving partial deriva- 
tives of .F13 to the first and second order, but it turns out 



Ati3 /l3,B Cp + /i3,c a + /i3,/3 b 



Ci 



(11) 



with 



ri3,x 



dLogJ^is 



a 
a 



dLogX 



and ^ 



(12) 



where tp is the characteristic duration of the pulse. For 
two given spectral bands and an assumed spectral shape, 
Eq.(ll) provides a linear estimate of the lag, which di- 
rectly shows how it is related to burst spectral evolution 
via the temporal derivatives oi Ep, a and /3. The "cur- 
vature parameter" \Ci\ depends on the pulse shape at 
maximum, large (resp. small) \Ci\ values corresponding 
to spiky (resp. broad) pulses. 

3. The lag luminosity relation 

Equation (11) gives the lag between BATSE bands 1 and 
3 if the values of Ep, a, (3, their time derivatives, and 
the pulse shape are known at maximum. It will become 
a LLR if these parameters can be related in some way to 
the luminosity. The Amati relation (Amati et al, 2002) 
provides such a link but, in its most studied version, it 
connects the isotropic energy in gamma-rays to the Ep 
value of the global spectrum. However, it has been sug- 
gested that a similar relation may exist between Ep and 
the luminosity. Yonetoku et al (2004) find, for example, a 
relation between the maximum luminosity and the global 
Ep, while Ghirlanda et al (2005) propose a relation be- 
tween the values of Ep and the luminosity both taken at 
pulse maximum 



i?p- 380(^g^^52erg.s-i 



0.43 



keV . 



(13) 



m) , (9) 



If a substantial fraction of bursts satisfy Eq.(13), it will, 
together with Eq.(ll), lead to a LLR that can be com- 
pared to the observational data. 

3.1. A (too) simple example 

In this section we limit the burst's spectral evolution to 
the variation of Ep only and write d = 6 = 0. Equation 
(11) then simply becomes 



At 13 /l3,£; Cp 



with fis^E = 



dLogJ^i3 



OLogEp 



ti,a,0 



(14) 



The function fis^E tends towards for the highest and 

lowest values of Ep, while it is at its maximum for 50 < 
Ep/{1 + z) < 100 keV. This is illustrated in Fig.l and 
can be easily understood, since for E <^ Ep (resp. E ^ 
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Fig. 1. The function /13.B given by Eq.(14) plotted as a function of the observed peak energy for a = —1 and 
P = —2.25; full line: linear scale; dashed line: logarithmic scale. 



Ep) the Band function B{x) behaves as x" (resp. x^). 
Therefore for high Ep values {Ep/{1 + z) :$> 300 keV) 

^13^ 501+0,201+" =^ ^""^ = ^ ^^^^ 

a similar result with /3 replacing a being obtained for 
low Ep/ (1 + 2) <C 20 keV. At intermediate values, 50 < 

Ep/{1 + z) < 100 kcV, and making the rough approxima- 
tion that B{x) oc a;" in band 1 while B{x) oc x^ in band 
3, we obtain 



3001+^ - 1001+^* / Ep 

^13 — 



50i+«-20i+« \l + z 
and therefore 



/: 



13. E 



■p. 



(16) 



(17) 



Finally it can be shown that the shape of B{x) at small 
X, B{x) = a;" [1 — (2 + a)x] leads to a power law behavior 
for /i3,_E at large Ep as seen in Fig.l. 

If the Amati-like relation (Eq.(13)) is satisfied, high 
luminosity bursts will have a large Ep and therefore a 
small lag, while lags will be comparatively large for bursts 
with Ep in the range 50(l+z) - 100 (1+z) keV, i.e. L ~ 
4 10^° (l-l-^;)^-^ erg.s at pulse maximum. Equation (14) 
also predicts that, for a given luminosity and spectral evo- 
lution, spiky bursts (large |Ci |) will have smaller lags than 
bursts with broad pulses (small |Ci|), in agreement with 



observations (Hakkila & Giblin 2006). The LLR obtained 
with this simple model is represented in Fig. 2 for differ- 
ent values of the ratio \ep/Ci \ and a typical burst redshift 
z = 1. A large (resp. small) ratio corresponds to a faster 
(resp. slower) spectral evolution or to a broader (resp. 
spikier) pulse. At high luminosity {L > 10^^ erg.s^^), the 
LLR has a power law behavior since 



dLogL _ / dLogL \ / dLogEp\ 
d Log At ~ \dLogEp) \d Log At J 



(18) 



the first factor being given by the Amati relation (Eq.l3), 
while the second results from the spectrum shape. From 
Eq.(14) we have 

d LogEp d LogEp 



d Log At dLogfi3^E 



(19) 



which is a constant at high Ep (and hence luminosity) val- 
ues (see Fig.l). However at luminosities smaller than 10^^ 
erg.s" ^ the model predicts that the power law behavior of 
the LLR should break down with the lag passing through 
a maximum and then decreasing. This clearly contradicts 
GRB 980425, which has both a very low luminosity and 
a large lag. But GRB 980425 docs not satisfy the Amati 
relation so that its departure from the LLR is not sur- 
prising. One should instead consider that this burst has a 
broad temporal profile, i.e. a small \Ci \ and an Ep of 138 
keV (Ghisellini et al 2006) corresponding to the maximum 
of fi3,E (see Fig.l) and therefore to a large expected lag. 
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Fig. 2. Lag-luminosity relation resulting from Eq.(14) which limits the spectral evolution to a decrease of Ep during 
pulse decay (a and /? being kept constant respectively equal to —1 and —2.25). The three lines cover an order of 
magnitude in |ep/Ci| from 0.1 (dashed line) to 0.3 (full line) and 1 (dotted line). A redshift value z = 1 has been 
assumed. 



3.2. A more complete study 

The spectral evolution of GRBs is, however, not limited to 
a decrease in Ep during pulse decay. A hard-to-soft evo- 
lution is also observed for the spectral indices a and /3. 
In some extreme cases a has been seen to decrease from 
a ~ 1.5 (a value a priori excluded by the synchrotron 
model) to about —0.5 in just a few seconds (Crider et 
al 1997). When the variation in the spectral indices is in- 
cluded in our linear analysis, it no longer predicts a vanish- 
ing lag for high or low Ep and L values since now /i3,a 7^ 
when Ep ^ 00 and /i3,/3 7^ when Ep ^ (for exam- 
ple limsp^oo /i3,a — -1.7 for a = — 1 and /? = —2.25). 
The lag then reaches a constant limiting value at low and 
high luminosities, where it apparently contradicts the ob- 
served LLR. This is shown in Fig. 3 where our calculated 
LLR has been plotted for different values of d and a fixed 
b = 0.1. Even a moderate variation in the spectral indices 
has a dramatic effect on the LLR and the global agree- 
ment with the Norris et al (2000) results that was found 
in the last section is now lost. 

4. Discussion 

These results clearly disagree with the observational data 
for the most luminous GRBs. Therefore if real bursts do 
satisfy the LLR proposed by Norris et al. (2000), a solution 



has to be found for the apparent discrepancy between our 
analysis and the observations: 

(i) A first option could be that in most cases the vari- 
ation in the spectral index a is small, at least around 
pulse maximum. However this does not seem to be the 
case for the bright events for which a detailed, time re- 
solved, spectral analysis has been possible (Precce et al 
2000). Moreover, the constraint on any variation in a ap- 
pears so severe (only the LLR with d = 0.01 in Fig. 3 is 
marginally compatible with the Norris et al. results) that 
it seems difficult to expect it will be satisfied by a large 
fraction of GRBs. 

(ii) A more interesting possibility would be that a re- 
lation may exist between the curvature parameter \C'i\ 
and the luminosity, bursts with spiky pulses being more 
luminous than bursts with broad pulses. This might be a 
different way to express the variability- luminosity relation 
proposed by Reichart et al (2001). We tried, for example, 
a simple linear expression of the form 

|(7i| = 2 + 0.2(^51-1) (20) 

where the resulting LLR with |ep/Ci| = 0.3 and d = b = 

0.1 is represented in Fig. 3. In spite of the variation in a, it 
now gives again very small lags at high luminosity because 
the pulses are then much spikier than at low luminosity. 
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Fig. 3. LLR with |ep/Ci| — 0.3, but now also including the variation in the low and high-energy spectral indices. The 
thin lines correspond respectively (from left to right) to d = 0.01, 0.03, 0.1, and 0.3 (and have all b = 0.1), while the 
thick full line represents the d = b = case. The dashed line is obtained with d = b = 0.1 and a varying curvature 
parameter given by Eq.(20). 



5. Conclusion 

We have performed a linear analysis of the time lag be- 
tween two spectral bands and have obtained a simple rela- 
tion (Eq.(ll)) which clarifies how the lag is related to the 
burst spectral evolution. When this relation is used in con- 
junction with the Amati relation it leads to a satisfactory 
LLR only if the spectral evolution of GRBs is limited to a 
decrease of Ep during pulse decay. If the variation in the 
spectral indices is also included, the lag does not decrease 
any longer to low values, even at very high burst lumino- 
sity. We have briefly discussed the possibility that short 
lags might be recovered if burst luminosity is correlated 
to the shape of the pulses, bursts with spiky pulses being 
more luminous than bursts with broad pulses. The avail- 
able sample of GRBs with both measured lags and known 
distance is still small but should increase with SWIFT. 
This will allow constraining tests of the results presented 
in this paper. 
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